Ground-State Degeneracy of Correlated Insulators with Edges 
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Using the topological flux insertion procedure, the ground-state degeneracy of an insulator on 
a periodic lattice with filling factor = p/q was found to be at least g-fold. Applying the same 
argument in a lattice with edges, we show that the degeneracy is modified by the additional edge 
density Ve associated with the open boundaries. In particular, we demonstrate that these edge 
corrections may even make an insulator with integer bulk filling degenerate. 



A fascinating aspect of quantum mechanics is the in- 
terplay between global topology and some physical prop- 
erties, that are naively believed to depend only on local 
quantities. Restrictions that originate from topological 
arguments [1] often lead to elegant and non-trivial pre- 
dictions of bulk properties in a non-perturbative fashion. 
One example is particle statistics dictating the possible 
ground states of a system. A stark demonstration of 
that is anyons with fractional statistics [2]. Another in- 
stance of this picture is the non-perturbative proof of 
Luttinger's theorem[3, 4]; it is revealed that the topol- 
ogy dictates the space enclosed within the Fermi surface 
of a Landau-Fermi liquid to be exactly the particle den- 
sity in the conduction band, disregarding the details of 
the microscopic Hamiltonian. Beyond their elegance and 
appeal, topological ground-state properties have been re- 
cently made the centerpiece, primarily due to Kitaev, of 
a novel kind of quantum computation[5]. 

Predictions from topological arguments are often sen- 
sitive to the geometries of the bulk. For instance, the 
insulating Z2-gauge spin liquid [6] is four- fold or two- fold 
degenerate depending on the bulk geometry is a torus 
or a cylinder respectively. Another well-known example 
is the V = 1/q fractional-quantum-Hall liquid [7] on the 
two-dimensional torus. A topological constraint due to 
these exotic excitations gives rise to the non-trivial q- 
fold degeneracy. Remarkably, by cutting the torus into 
a closed strip with two open edges, the ground-state de- 
generacy in the thermodynamic limit shoots up from q to 
infinity, with gapless excitations described by the chiral 
Luttinger liquid on either side of the Hall bar. It is then 
natural to expect that the topological arguments might 
deliver rather different messages when edges are present. 
Indeed, the quantum Hall system is not the only example 
of juicy edge physics. Other examples are Andreev bound 
states appearing on edges of superconducting lattices [8], 
spin- 1/2 excitations of the Haldane phase located at the 
edge of S* = 1 Heisenberg chain [9], and the ferromagnetic 
moment on a zigzag tip of a carbon nanotubc[10]. 

In this Letter, we revisit a topological constraint 
in insulators - the correspondence between the fill- 
ing factor and the ground state degeneracy. In 



Ref. [11], Oshikawa extends the Lieb-Schultz-Mattis 
argument [12] and Laughlin's treatment of the quantized 
Hall conductance [13], to show that the ground state de- 
generacy of an insulator with filling fraction v ~ p/q 
(where q and p are coprimes) is at least g-fold. This 
is established using the flux-insertion method applied to 
the system on the periodic and finite lattice. By taking 
appropriate thermodynamic limit, the ground state de- 
generacy is related to the intensive particle filling factor 
in the bulk. However, in the presence of open boundaries, 
the particle density is no longer uniform near the edges. 
The physics of the edge is captured through appropriate 
edge filling factors that are introduced in addition to the 
bulk density. In this Letter, we will show that the degen- 
eracy of the ground state is determined by both the hulk 
and edge filling factors in the presence of edges. 

For instance, a gapped insulator on a periodic lat- 
tice, made of spinless particles with bulk filling factor 
i'b = 1, may have no degeneracy. However, if the lat- 
tice has open boundary conditions along one direction, 
the particle density near the edge may differ from the 
bulk filling factor far from the edge. This difference is 
captured by the edge density v^, which arises from par- 
ticles localized near the edge, as well as changes of the 
extended wave functions in the bulk. The edge state may 
break translational invariance and gives rise to ground- 
state degeneracy despite of the bulk state having perfect 
translational symmetry parallel to the edge. For exam- 
ple, for Ve — 1/2, the edge may form a CDW with period 
two. Such a state is doubly degenerate, although the bulk 
insulator is featureless. 

Let us now review the flux insertion procedure in the 
context of a system with open boundaries in one direc- 
tion (we follow Ref. [11] closely). Consider a lattice in 
d-dimensional space with open boundaries in one direc- 
tion, 1 < cci < Li, and periodic in the remaining d-\ 
dimensions with finite length Li, i — 2, . . . ,d. The insu- 
lating ground state |$o) contains N particles that reside 
on the lattice sites. Making use of the periodicity of the 
lattice, say in the Xd direction, the ground state can be 
chosen to be an eigenstate of the momentum operator 
Pd, carrying momentum in the Xd direction. 



2 



We now introduce a fictitious vector potential A that 
couples with "unit charge" to all particles in the system. 
Protected by the gap, one can adiabatically insert a flux 
quanta ^ = he through the ring in the .x^i-dircction and 
map |\E'o) to of the same energy. Because the flux 
insertion can be achieved by the constant vector poten- 
tial Ad = ^/Ld in the direction of .t^, which commutes 
with Pd, the momentum remains constant in the whole 
adiabatic procedure, Pd\%) = -P°|*o)- 

But the flux insertion changes the Hamiltonian from 
7i($ = 0) to a different topological sector 7i{^ = he). 
Before being able to compare to |4'o), we need to 
restore the Hamiltonian to the same topological sector 
= 0). This is done using the following large unitary 
gauge transformation: 

U = exp ^ Xdhs^ . (1) 

Now Un{hc)U-^ = n{0), and therefore C/l^-^) is an 
eigenket of the original Hamiltonian. Since our system is 
a gapped insulator, and the flux insertion was adiabatic, 
we expect that UI^^/'q) has the same energy as |\E'o) (with 
n{0) in mind). 

The momentum of the newfound ground state can 
be evaluated straightforwardly, P^C/l^'o) ~ {UPd + 
[Pd,U])\%) = {P2 + '2nNlLd)U\%). So the momen- 
tum of U\%) is P'd = PS + 2TrN/Ld. If N and Ld are 
mutually prime, then U\^q) and |^'o) are two degener- 
ate ground states. To relate the ground state degeneracy 
to the fllling factor, one need to take appropriate ther- 
modynamic limit. In a fully periodic system, Oshikawa's 
argument follows by saying that 

AP = - = 27r^ = 27ri/Q, (2) 

where i/ = N/V = p/q is the filling factor of the lat- 
tice. The volume of the system is ^ = ni=i ^^'^ 
transverse size at each particular Xi is Ct = V/Li. The 
thermodynamic limit is taken by fixing v (so that p and q 
are well defined without any correction) and choosing Cd 
to be coprime with q. Thus, by repeating the flux inser- 
tion, we can generate g-fold distinct ground states. If we 
assume that the ground state degeneracy Dg is a robust 
quantity, disregarding the details of the thermodynamic 
limit, the gauge argument leads to the constraint Dg > q. 

Returning to the non-periodic system introduced 
above, we pick up the discussion from Eq. (2). Two 
different kinds of non- uniformity occur near the edges. 
First, the edges give rise to new states of particles lo- 
calized near the edges; we deflne the edge-state filling 
factor Ue = Ne/Ci, where is the number of edge par- 
ticles. Second, bulk states (forming a continuum) that 
propagate throughout the lattice may be repelled from 
the edges. This repulsion creates a charge depletion near 
the edge; we deflne = N^/Ci, where is the total 



charge depleted near the edge. Both effects give rise to 
non-trivial corrections to the gauge argument. 

Loosely speaking, the presence of edges divides the to- 
tal number of particles into two groups: bulk and edge, 
labeled by Ng and respectively. The precise distinc- 
tion between bulk and edge states would be explained 
in later paragraphs. In light of this division, we need to 
restate Eq. (2) as 

AP = 277 ^%+^" = 27r^ + 2nneCu, (3) 

with Cid = V/{L\Ld). Before we can conclude anything 
about the ground-state degeneracy, we also need to ac- 
count for the repulsion of bulk states from the edge. 

As mentioned above, the edges may repel the bulk 
states, creating a charge depletion near the edges with 
density = N^/Ci. Now, in order to achieve the ther- 
modynamic limit of an insulator with bulk filling i/b said 
edge filling Ug, we need to consider a sequence of finite 
lattices with a total number of particles which is 

No = Ng + N^ = {ugV - N^) + n^Cx 

= UbV + {ne-n^)Cx (4) 

Note that tIa = P^/q^ is not necessarily an integer. To 
constrast, in a periodic system the thermodynamic limit 
would simply be taken by considering a sequence of lat- 
tices with Nq = VbV particles. 

Including both types of edge effects, the momentum 
difference between and its sibling U\^'q) is 

AP = 2T:VgCd + 2ttveCu. (5) 

with Vg — Pslqs and — — ~ Pm/qE- From the 
revised gauge argument in Eq. (5), we obtain the main 
result of this paper: the degeneracy of an insulating state 
in the presence of edges is given by 

Dg>LCD{qn,q^)>qg, (6) 

where LCD denotes least common denominator. 

Before discussing possible extensions and generaliza- 
tions of the result in Eq. (6), we resolve some subtitles 
in the derivation above. A central issue is the difficulty 
of distinguishing edge and bulk particles as in Eq. (3) in 
a strongly correlated system, where the concept of single 
particle states is rather vague. At first sight, this dis- 
tinction may seem impossible to carry out in practice, 
and might be mistaken to be a theorist's whim; there is, 
however, a precise way of counting the number of edge 
and bulk states. This method utilizes a spectral fimc- 
tion decomposition of the system's correlation function, 
which in principle we know with utmost precision. Fol- 
lowing Lehman decomposition of the single-particle spec- 
tral function: 

A{x,u:) = Y,^{u - E^)\{m\i,^{x)\^o)\'' 

m 

± ^<5(a; + i;„)|(m|V(x)|*o)|'. (7) 

m 
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Here ^pix) is the annihilation operator in the Schrodinger 
picture at lattice site x. The ± distinguishes between 
fermions and bosons. Without loss of generality, let us 
concentrate on fermions. Note that |to) are many-body 
excited states that show one particle excitations. The 
excitations correspond to one-particle states of free par- 
ticles. In the following we simply refer to |m) as states. 

Next, consider the a; < (hole) part of the spectral 
function A{x,u)). In a finite system, it is given as a dis- 
crete sum of 5-functions multiplied by matrix elements, 
and it obeys the sum rule 



j A{x,i^)du} = N, (8) 

^ — oo 

where N is the total number of particles in the lattice. If 
we take the limit Li — > oo while keeping the bulk density 
fixed, we expect that some of the (5-functions will merge 
into continuum, whereas the rest will remain sharp and 
seperated. The first group is the excited hulk states, and 
the second group is the excited edge states. We stress 
that |m) are strongly correlated many-body states that 
show one particle excitations AQ = ±1 relative to the 
ground state. More precisely, bulk states Ims) scale as 

|(m3|V^(f)|*o)r~o(;^) (9) 

for any x, as Li is scaled to infinity. Unlike bulk states, 
the edge states, Ime), tend to an a?-dependent constant 
as L\ — > oo: 

\{m^\^l^\x)\^o)? 0{l). (10) 

Thus we can break the spectral function into a bulk {Ab) 
and edge {Ae) parts: 

± Y 5{^ + EmB/Mmn/.\i^{S)\^^)?. (11) 

mB/E 

The definition of Nb and Ne is hence 



NBiE) = Y j di^AB^E){x,uj). (12) 

^ — oo 

There are two subtleties in the above procedure. First, 
it is possible that an edge state hybridizes with a bulk 
state. Definition along the lines of Eqs. (11) and (12) 
specifies that such a state is an edge state; even though 
such a state is partially delocalized, it is associated with 
the edge due to the localized weight. The number of 
such states will scale with Li for i > 2, but not with 
Li. Another subtlety occurs if there is an accidental 



degeneracy between a bulk and edge states. In this case, 
in a finite system, it may be impossible to distinguish 
between the edge and bulk states and two hybridized 
states. This accidcntel degeneracy will surely be lifted 
by a different choice of system dimensions, and therefore 
needs not bother us. 

Now that the derivation is complete, we would like 
to elaborate on some extensions of the edge argument. 
The first question we address is the effect of bulk-state 
depletion near the edge. Note that if is an inte- 
ger, the bound on the ground-state degeneracy, Dg > 
LCD{qB,qe), only depends on Vb and rig = Pe/Qe- In this 
case, the depletion of bulk states does not play any role 
in determining ground-state degeneracy. In fact, for fea- 
tureless insulators with low-energy excitations described 
by Fermi liquid theory, we can show that is an integer. 

Start with the ground state |0) of the featureless in- 
sulator on the periodic lattice (without edges). The in- 
sulator with open edge can be viewed as perturbing the 
original ground state |0) with hole excitations near the 
boundary. The translational invariance ensures that we 
can construct local quasi- hole states \r) = tp{r)\0), with 
the same spatial profile but located at different Ci lat- 
tice sites on the edge. In general, these states, located 
at different lattice sites, would have non-vanishing over- 
laps. If one assumes the Fermi-liquid picture is at work 
here, these local (Wannier) orbitals would form a con- 
tinuous band (with dispersion needs to be determined 
self-consistently). To ensure the new ground state (in 
the presence of open boundaries) is also insulating, these 
quasi-particles must form a band insulator at the edge, 
i.e. the number of holes in one unit cell, = N^/C\, is 
an integer. 

The argument also applies to a simple but exotic sit- 
uation when the low-energy excitations carry fractional 
"charges" (not necessarily the ordinary electric charge), 
but are well described by Fermi-liquid like theory. Fol- 
lowing similar argument, Ua is fractional, reflecting the 
breaking up of ordinary particles. While the argument 
does not hold when low-energy excitations are no longer 
described by Fermi liquid theory, in the above case, 
serves as a useful precursor for identifying low-energy ex- 
citations in featureless insulators. 

Another extension of the theorem is when the edges 
are far apart. It is then natural to divide the edge cor- 
rections. He = and = n^^, associating with 
each edges. A simple example helps to demonstrate this 
point and the statements above; consider a segment of 
carbon nanotube with circumference Ly (in y direction) 
and zigzag edge at a; = and x = L^. Choose the tight- 
binding hopping to be in the strong anisotropic regime, 
\th\ > 2\t\, where th and t denote the hopping amplitudes 
along horizontal and (vertical) zigzag bonds. The inter- 
actions between particles are assumed to be weak, with 
a typical energy scale V . 

Ignoring the weak interaction momentarily, the open 
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boundaries at x = and x = L^, give rise to single- 
particle edge states at -E = with momentum-dependent 
localization length ^{ky) = In \ {2t cos ky)/th\- Each edge 
has Ly distinct states, corresponding to each quantized 
ky. Since the bare hopping Hamiltonian is quadratic, the 
total number of states has to equal the number of parti- 
cles in the lattice, or twice the number of unit cells (each 
carbon-nanotube unit cell contains two sites). Hence, 
the Ly edge states must come at the expense of Ly bulk 
states. This leads to = 1 at natural filling i^g = I, 
and thus there is no effect on the degeneracy due to the 
bulk states repulsion at the edge. 

Nevertheless, the edge-state filling factor, = Uei + 
ne2, rnay be tuned to arbitrary fractions by removing or 
adding a small amount of electrons (proportional to the 
length Ly) which do not affect Vg- Because the interac- 
tion V is weak (compared with the bulk gap A^), the 
non-trivial mixing of single-particle state occurs within 
edge states and lifts the exact degeneracy at £" = 0. The 
ground state is non-trivial and certainly depends on the 
specific form of interactions. But the revised gauge ar- 
gument with the momentum shift in Eq. (5), leads to at 
least Qe-iold degeneracy of the ground state even though 
the bulk filling is an integer. 

When considering particular physical systems, as in 
the above example, it is possible that the ground-state 
degeneracy will be larger than our rigorous result in 
Eq. (6). If the two edges are known by other means 
to be independent (as is probably the case in many 
large systems) and are also independent of the bulk 
state, then the degeneracy should be at least the product 
gel • Qe2 ■ LC D{q^,qB). This degeneracy would become the 
number of low- lying states, if there is a weak interaction 
between the two edges. Other cases are also possible, 
most notably, independence of the edge states from the 
bulk, but not from each other. In this case the degen- 
eracy is LCD{qei. ■qe2) • LC'D{q^,qB). We add that for 
inconirnensurato edge filling in the thermodynamic limit 
i — > oo, it is very likely that the low-energy physics of the 
system is described by a gapless liquid on the edge, liv- 
ing inside the gapped bulk spectrum (similar to the chiral 
edge states in quantum Hall liquid). Since these edge ex- 
citations sometimes are the only low-energy excitations, 
appropriate treatment of them is crucially important. 

An interesting situation that may also be quite com- 
mon is that the bulk and edge conspire to produce a uni- 
form particle density even near the edge. This happens 
if the bulk states repelled from the edge are exactly com- 
pensated by particles trapped at the edge: rig = n^. Our 
revised gauge argument is then reduced to Oshikawa's 
original argument with only the bulk filling determining 
the degeneracy. 

The revised gauge argument can be generalized to the 
case of more than one non-periodic direction. We require 



that at least one of the directions of the lattice is peri- 
odic to derive a generalized gauge argument, although 
this restriction may not be necessary in a physical sys- 
tem. We consider Xd as periodic. If all other sides of 
the attice terminate at = and Xi = L, there may 
be edge (surface) and wedge states, with filling numbers 
that scale as ^ i", where a = 1, 2, 1. One then 
needs more filling factors Ua to specify the filling of the 
system. Also, one needs the depletion parameters 
where a = 1, . . . ,d specifies the repelled 'bulk' state, and 
b = l,...,a — 1 is the number of restricted dimensions 
for the bulk states. To clarify this statement consider 
the example of a cube that is made periodic in the z 
direction. The bulk is three dimensional, and the deple- 
tion away from the surfaces {x = 0, x = L, y = and 
y = L) amounts to n^j ■ L^ particles. Also, the wedges 
at X = 0,y = 0, X = 0,y = L etc. may repel the bulk 
states with a volume of n^^ ■ L. Similarly, edge (surface) 
states on the faces of the cube may be repelled from the 
wedges with corrections n]^^ • L. 

In summary, using topological flux insertion procedure, 
we rederive the connection between the ground-state de- 
generacy in an insulator with its filling factors i/g in the 
bulk and i^e at the edge. As expected, the presence 
of open edges produces rich physical phenomena at the 
boundaries. In addition, the response of the bulk states 
to the edge may lend a different perspective into the cor- 
related bulk physics. 
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